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SUMMARY 

Equations for a lumped parameter mathematical model of a subsonic wind ‘ 

tunnel circuit are presented. The equation state variables are Internal energy. / 

density, and mass flow rate. The circuit model Is structured to allow for 
integration and analysis of tunnel subsystem models which provide functions 
such as control of altitude pressure and temperature. Thus the model provides 
a useful tool for Investigating the transient behavior of the tunnel and con- 
trol requirements. The model was applied to the proposed NASA Lewis Altitude 
Wind Tunnel (AWT) circuit and Included transfer function representations of 
the tunnel supply/exhaust air and refrigeration subsystems. Both steady-state 
and frequency response data are presented for the circuit model Indicating the 
t^pe of results and accuracy that can be expected from the model. Transient 
daca for closed loop control of the tunnel and Its subsystems are also 
presented, demonstrating the model's use as a controls analysis tool. 


INTRODUCTION 

As a rcraft propulsion technology continues to advance to meet more demand 
Ing mlss'on requirements, the need for Improved and expanded testing facilities 
will Increase. For example, some future aircraft will have more highly Inte- 
grated propulsion and airframe systems. The testing facilities for these 
future aircraft will need to have test sections that can handle full size 
models. Other requirements will be to concurrently simulate pressure and tem- 
perature to accurately reflect altitude conditions, to have air speeds 
approaching Mach 1.0, to provide good flow quality and anecholc properties In 
the test section, to have operating propulsion systems In the test sections, 
and to adequately simulate adverse weather conditions (particularly heavy rain 
and Icing). To Insure a technically sound design for these new facilities, 
much physical modelling (large and small scale) and analytical modelling and 
computer simulation will be necessary. This report describes a lumped param- 
eter mathematical model for simulating a subsonic wind tunnel circuit. The 
model was developed at Lewis. The model and Its resulting simulations (digital 
and/or analog) can be used as a test bed to determine system dynamic response 
characteristics, to assess control design concepts, to specify operating sce- 
narios, and to establish an analytical database. Additionally, various sub- 
systems can be Integrated Into this model to further expand the above analyses. 

The set of nonlinear differential equations for a five-lump circuit model 
are presented, together with simplified transfer functions that can be used to 
represent tunnel subsystems (e.g., refrigeration and air supply/exhaust). 
Steady-state and dynamic (l.e., frequency response) results are presented for 
a specific wind tunnel design to demonstrate the capabilities of the model. 
Transient results are also presented to demonstrate the use of the model for 
controls analyses. 
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ALTITUDE WINO TUNNEL APPLICATION 


The lumped parameter model has been applied to the analysis of a proposed 
rehabilitation of the AWT at Lewis. The proposed rehabilitation would provide 
the aforementioned aerodynamic propulsion test capabilities. Reference 1 gives 
a more complete description of the tunnel and Its proposed capabilities. The 
AWT wind tunnel circuit with Its major subsystems Is shown In figure 1. The 
fan drive system circulates air In the tunnel and Is the means for controlling 
test-section Hach number. The function of the supply and exhaust air systems 
Is to establish the desired altitude pressure condition. Air Is exhausted from 
the tunnel to lower pressure (raise altitude) or Injected Into the tunnel to 
raise pressure (lower altitude). The refrigeration system Is used to remove 
the fan heat of compression and to maintain the desired altitude temperature, 
finally there Is a plenum evacuation system (PES) that Is used to bleed the 
boundary layer out of the test section to minimise tunnel wall Interference 
effects. Following paragraphs will describe the circuit model along with the 
models of the fan drive, supply/exhaust air, and refrigeration subsystems. 


TUNNEL CIRCUIT MOOFL 

The aerodynamics of a wind tunnel circuit Is correctly represented as a 
three-dimensional distributed process. That Is. the fluid states are functions 
of both time and space. However, the solution of the associated partial dif- 
ferential equations Is computationally difficult and time consuming. In 
dynamics and controls applications, the need for rapid solution and reasonable 
dynamic fidelity override the need for extreme steady-state accuracy, for 
example, real-time solutions are often desired to provide a realistic means 
for checkout of control software. It was -In this light that a one-dimensional 
lumped parameter approach was selected to represent the Altitude Wind Tunnel 
(AWT) circuit. The lumped parameter approach approximates the distributed 
process by dividing the circuit Into a number of lumps, each covering a spec- 
ific length of the tunnel. The fluid states are then calculated In each lump. 
Using more lumps Increases the model's dynamic fidelity but also Increases Its 
complexity. The number and location of lumps are also chosen to facilitate 
Incorporation of the tunnel subsystems (e.g., supply/exhaust air and 
refrigeration) . 

for the AWT application, the wind tunnel circuit was divided Into five 
lumps as shown In figure 2. The five lumps provided a good compromise between 
capturing the Important dynamic characteristics of the tunnel and the model 
complexity (number of states). The station numbers shown In boxes In figure 2 
correspond to stations In a contractor-supplied computer program that makes 
detailed steady-state performance calculations (hereafter referred to as the 
steady-state performance deck). 

For each lump, ordinary differential equations were written to compute the 
time rate of change of mass flow rate (over the length of the lump) and the 
time rate of change of density and energy storage (In each of two half-lump 
volumes). These equations represent conservation of momentum, mass, and 
energy, respectively. The equations have no spatial dependence as would be 
the case for a distributed parameter model. Since the lumps are contiguous, 
adjacent half-lump volumes can be combined Into one flow volume with the time 
rate of change of mass and energy computed In each of these flow volumes. Note 
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that this combining of flow volumes effectively translates them spatially from 
the lumps (see fig. 3). 

The circuit model formulation Includes a fan. The purpose of the fan Is to 
circulate air In the tunnel, providing a pressure rise to compensate for the 
pressure losses around the tunnel circuit. These pressure losses are due to 
corner turning vanes, screens, diffusers, the cooler, etc. Static pressures 
are calculated from auxiliary equations In the volumes Immediately preceding 
and following the fan. From these pressures, a total pressure ratio Is cal- 
culated. Using this pressure ratio and a specified fan speed, the airflow 
through the fan Is obtained using a table lookup. Because the fan Is Included 
In the model, the two half-lump volumes on either side of the fan cannot be 
combined. Hence. If there are n lumps, there will be n ♦ 1 volumes In the 
model. Therefore, there are six flow volumes In this five lump model. 


The general form of the state equations Is given below. The state vari- 
ables In this lumped parameter model are Internal energy, density, and mass 
flow rate. There are 17 states in the five lump circuit model (no subsystems 
or controls). All symbols are defined In appendix A. 
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Equations (1) and (2) are one-dimensional dynamic forms of the conservation 
of mass and energy equations. The conservation of momentum equation (eq. (3)) 
assumes a constant area duct for each lump. These equations allow for the 
addl t lon/remova 1 of mass and heat Into/from the tunnel circuit. It Is also 
assumed, In equation (C), that pressure loss Is proportional to pV?. Included 
In equation (3) Is a term fj. This term Is referred to. In this report, as 
a frictional loss factor. It Is used to compensate for all pressure losses 
around the tunnel. At each steady-state operating condition of the tunnel, 
there Is a set of frictional loss factors which balance the momentum equations 
(l.e., cause m^ » 0). Each of these frictional loss factors could vary with 
the tunnel operating condition. Hence, It may be necessary to Include these 
terms In some functional form If the model Is used to Investigate large per- 
turbations from the operating condition. Some source of steady-state data 
(analytical or experimental) Is required to allow calculation of the "exact" 
frictional loss factors. 


For the AWT application, analytical data were available from the steady- 
state performance deck. The data Indicated that the term 

(? 1 ( L/0) ) (p^v^/ 2) g c J was nearly linear with (1/2) p,v^ for 17 conditions 
around the tunnel's operating envelope. Therefore, constant values for the f)’s 


3 



were determined from the slope of a least squares, first-order curve fit of the 
data. Results will be presented using both the ‘exact* and least square values 
of 7^ . Table I gives the values for both sets of frictional loss factors. 

The tunnel volume can be thought of as being composed of a series of frus- 
tums of cones. In order to provide a more accurate approximation of the flow 
volumes, each lump was subdivided, as described In appendix B. Since constant 
diameter lumps are assumed for each lump, changes In tunnel diameter over the 
length of. a lump were accounted for by precalculating volumes, lengths, 
weighted average areas, and weighted average diameters. Table II contains the 
values used for these areas, lengths, diameters, volumes, etc. In the AWT 
model. Detail on these calculations Is also provided In appendix B. 

Auxiliary equations are used to calculate the pressures, temperatures, and 
Mach numbers of Interest around the tunnel. The general form of these equa- 
tions Is as follows 
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Following Is the general form of equations for the fan. 
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Constant adiabatic efficiency Is assumed for the fan. The data used In the 
calculation of fan corrected airflow (eq. (10)) Is taken from a fan map such as 
the one shown In appendix C, figure i. Map data for the AWT fan are given In 
appendix C, table I. The complete set of equations for the circuit, Including 
the fan equations. Is given In appendix C. 


SUBSYSTEM MOOELS 

A schematic of the tunnel circuit and the major subsystems (l.e., fan drive, 
supply/exhaust air, and refrigeration subsystems) represented In the AWT appli- 
cation Is given In figure 4. The proposed control architecture which maintains 
the desired test section Mach number, temperature, and pressure Is also shown. 

In actuality, the subsystems are extremely complex processes that provide 
Inputs (e.g., flows) to the circuit and In general would receive outputs (e.g., 
measurements) from the circuit, as represented by the solid and dashed arrows 
In the figure. For the AWT application, the subsystems are simply represented 
by transfer functions (to be described) to demonstrate the use of the circuit 
model. However, the subsystems could be represented In much greater detail, If 
desired. Proport Iona 1 -plus- Integra 1 control logic was chosen for this study to 
achieve zero steady- state error In the tunnel control parameters. As can be 
seen from figure 4, the controls also receive Inputs from the tunnel circuit 
and supply actuator commands to the appropriate subsystems. 
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Fan Drive Subsystem Model 


A block diagram of the fan drive system and the Mach number controller Is 
given In figure 5. The Mach number controller operates on the error In the 
Mach number - the difference between the set point and the measured values. 

In actuality, Mach number would be calculated from appropriate pressure measure 
ments and a calibration curve. The output of the Mach number controller Is a 
speed set point signal which is fed to the speed control loop. A proportional 
speed controller Is used with Its output being a motor torque demand. The 
motor torque response Is assumed to be a first order lag which Is limited. The 
calculated fan torque Is subtracted from the motor torque and the torque dif- 
ferential (l.e., fan acceleration) Is Integrated to give fan speed. The values 
that were assumed for the various parameters In the AWT system are given In 
figure 5. 


Supply/Exhaust Air Subsystem Models 

A block diagram of the tunnel pressure control and the supply/exhaust air 
subsystems Is given In figure 6. Although test section total pressure Is shown 
as the controlled variable. It Is more likely that a pressure upstream In the 
bellmouth would be used In conjunction with a calibration curve. The output of 
the pressure controller Is a supply/exhaust air flow setpoint. A positive set- 
point Indicates that an Increase In pressure Is desired, resulting In a supply 
air demand. A negative setpoint Indicates that a decrease In pressure Is 
desired, resulting In an exhaust air demand. The flow-split logic passes the 
flow setpoint to the appropriate subsystem. Closed-loop flow control of the 
supply and exhaust air subsystem are simply represented by second order trans- 
fer functions. The assumed values for the second order natural frequency and 
damping were based on past experience with similar large facilities. It was 
assumed that the circuit has no effect on the operation of the supply/exhaust 
air subsystems. This Is a good assumption for the supply air subsystem becaut : 
It Is expected that the supply nozzles will be choked. The assumption Is a 
simplifying one In the case of the exhaust subsystem. A change In tunnel pres 
sure would affect the amount of air exhausted unless the exhaust subsystem 
control valve Is adjusted. Assumed values for the various subsystem constant 
parameters are given In figure 6. 


Refrigeration Subsystem Model 

A block diagram of the tunnel temperature control and the refrigeration 
subsystem Is given In figure 7. The temperature control variable Is shown as 
test section total temperature. But, as with pressure, the temperature would 
probably be measured at an upstream location, such as In the bellmouth. The 
refrigeration subsystem Is a much simplified version of a more detailed model 
(unpublished). A brief explanation of the operation of the refrigeration 
system and Its representation In figure 7 Is as follows. Liquid freon from a 
storage tank known as the suction trap Is pumped through the heat exchanger 
colls located In the tunnel. The rate of heat extracted from the tunnel air, 
Q HX , *s proportional to the difference between the temperature of the colls, 

T CO I L * ancl the average temperature of the air flowing across the colls, 
T AIR,AV6- The heated freon Is returned to the suction trap as liquid and vapor 



The temperature of the liquid freon In the suction trap Is maintained by pump- 
ing the freon vapor out of the suction trap, through a cooling tower, and 
returning It as liquid to the suction trap. The closed loop control of the 
liquid freon temperature, T^sj. Is represented In figure 7 by a simple first 
order lag transfer function. The output of the tunnel temperature controller 
Is the temperature setpoint for the liquid freon In the suction trap. The rate 
of change of the temperature of the heat exchanger colls depends on the dif- 
ference between the heat load on the tunnel, $hx* and the amount of heat 
extracted from the freon, Qfreon. shown In figure 7. The values of the 
parameters used In the AWT control and refrigeration subsystem are given In 
figure 7. Note that 3 of the parameters are a function of the operating 
condition. 


DIGITAL COMPUTER IMPLEMENTATION 

The model was programmed In Fortran and was executed on the Lewis 
IBM 370/3033 computer system. A second order Runge Kutta Integration scheme 
was Implemented, using a time step of 0.01 sec. The simulation consists of an 
Initial condition program, a transient program, and a plotting program. The 
transient program package forms the main part of the simulation. It contains 
the numerical Integration and store/print subroutines. !n this simulation, the 
Integration routine calls the subroutine containing the circuit loop and sim- 
plified subsystems equations. Because of this, the Integration method used can 
be easily changed. Data Is passed between the three programs In this simula- 
tion package through COMMON blocks. These COMMON blocks overlay when the next 
program In the package Is loaded Into core. 

The Initial conditions on the state variables are set In the Initial con- 
dition program. Other tunnel variables and parameters are also Initialized In 
this program. Additionally, In the Initial condition program, setpoint values 
for the test section operating condition can be entered. Changes to gains, 
time constants, etc., of the subsystem controls can be entered In this program. 
This Input Is through namelists. In the transient program, a namelist Is used 
to enter transient duration, time step, number of plot points stored, etc. The 
plot program takes the stored plot points and generates the desired plots. If 
a steady-state listing of selected engine variables Is desired, another program 
Is used to obtain this. Other packages such as a frequency response program 
can be run using the plot data. 


HYBRID COMPUTER IMPLEMENTATION 

In order to more readily evaluate tunnel controls, the tunnel model was 
also Implemented on a hybrid computer. The hybrid facility at the Lewis con- 
sists of two PACER 600 hybrid systems, manufactured by Electronics Associates, 
Inc. Each system consists of a digital minicomputer, two analog computers and 
an Interface, which provides for control and setup of the analog computers by 
the digital computer and data exchange between the analog and digital com- 
puters. A multivariable function generation (MVFG) system Is also available 
for use by either or both hybrid systems. It was originally planned to Imple- 
ment as much of the tunnel circuit model and the faster subsystem models, such 
as the process air system, as possible on the analog computers. Slower sub- 
systems, like the refrigeration system would be programmed on one of the digital 



computers. The second digital computer was to be used for the Implementation 
of digital control algorithms and for data collection and display. Both digital 
computers would be also used to set up and check out the analog computed models. 
Multivariable functions, such as fan and compressor maps, would be programmed 
on the MFVG system. Since the tunnel and subsystem models were quite non- 
linear, attempts were made from the beginning to minimize the number of actual 
multiplications and divisions performed on the analog computers. It was antic- 
ipated that the limited number of multipliers on the analog computers (24 or 
30/console) would dictate how much of the overall model, tunnel and su systems, 
would have to be programmed on the digital computers. Consequently, some of 
the model equations were rewritten and/or combined and a few approximations 
made to reduce -the number of multipliers required. 


The first approximations were made In the calculation of the tunnel total 
pressures at the various stations. The maximum Mach number In the test section 
was assumed to never exceed 0.9, l.e.. MTS <0.9. The approximation used for 
the test section total pressure was 
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with a maximum error of 0.55 percent. The Mach number In the remaining model 
stations was assumed to never exceed 0.4. The approximation used for the total 
pressure In the remaining stations was 

P T = P (1 ♦ 0.2 M2) 7 /? a P ( 1 ♦ 0.7 M?) (19) 

with a maximum error of 0.41 percent. 


Another approximation was made In the fan total temperature equation. The 
expression. 
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was fit, over the range of the pressure ratio, by a curve which was then 
approximated by 15 straight line segments and programmed on a digitally con- 
trolled function generator. The accuracy of this function was within the 
accuracy of the function generation hardware, l.e., 0.5 percent. 


To save several analog multipliers, the equations for static and total 
pressure, static and total temperature, and Mach number for each station were 
combined such that the actual equations programmed were the following 
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Static temperature was not computed at any station and the Mach number was 
calculated only In the test section. The remaining tunnel equations were 
programmed as written. 

A programmatic decision was made to use simplified versions of the sub- 
systems for tunnel control studies. As a result. It was possible to mechanize 
the subsystem models, and their respective controls, entirely on the analog 
portion of the hybrid system. Using these simplified subsystems and with the 
previously mentioned component-saving assumptions In the tunnel model, it was 
possible to fit all models on two analog computers. The digital portion of 
the hybrid system was used only for setup and checkout of the analog computers 
and for data collection and display. 


Once the tunnel model was Implemented on the hybrid system, It exhibited a 
tendency to drift slowly when run open-loop (without controls). The only way 
to prevent this was to close the process air control loop. It Is felt that 
this drift was attributable, at least In part, to the Integral nature of the 
closed-tunnel flow process. This drift was not considered to be a serious 
problem, since the actual tunnel would behave in the same manner If it were not 
perfectly balanced ( 1 . e . , changes In temperatures and pressures would occur If 
no control was maintained). 


MODEL VERIFICATION 

For the simulation to be used with confidence as a controls test bed. It 
had to adequately represent the actual tunnel behavior over a range of condi- 
tions both In steady-state and dynamically. To validate the Integrated circuit 
simulation (with subsystems and controls) several means were used. Since no 
experimental dynamic or transient response data were available from the AWT 
circuit. It was necessary to find another means of verifying the lumped param- 
eter approach. Since the wind tunnel Is actually a distributed parameter 
process, one means of verification would be to compare results from the lumped 
parameter model to results from a distributed parameter model. It was felt 
that good agreement would give credence to the lumped parameter model. 


Steady State Results 

Steady state data were obtained from both the digital and the hybrid com- 
puter simulations. Six operating points were selected to cover the operating 
envelope of the tunnel. Figure 8 shows the operating envelope for the tunnel 
with the six selected operating points circled (l.e., points 3, 4, 8, 9, 13, 
and 16). The test section conditions corresponding to all points or, the enve- 
lope are listed In table III. At each of the selected points, two sets of 
steady-state data were obtained. One set had the frictional loss factors set 
to the values that were calculated at the design point, point 13 (Mach 
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number - 0.8 and altitude - 32 000 ft). The second set had the frictional loss 
factors set to the values determined using the least squares method described 
earlier. In all cases, the temperature, pressure, and Mach number control 
loops were closed. Table IV shows the percentage error between these steady- 
state data and data obtained from a contractor-supplied steady-state perform- 
ance deck. The fan variables that were compared are speed, pressure ratio, 
rate of heat addition, and the torque required. Also Included Is an overall 
comparison o' the total pressures, total temperatures, and mass flow rates In 
all tunnel volumes (lumps). As can be seen In table IV, the percentage errors 
In m, Pj, and Tj are all less than 4 percent and were considered satis- 
factory for control study applications. 

Errors In the fan variables are seen to be larger, particularly for 
and T FAN- The reason Is that these variables are proportional to the change 
In total temperature, ATj, across the fan. Since the fan Inlet temperature Is 
held constant by the controls, small errors In fan. exit temperature can result 
In large errors In the aTj, and consequently In Qp^ N and tc^ n . For 
example, at point 13 conditions, a 1 percent Increase In the fan exit temper- 
ature results In a 25 percent Increase In the alj. Table IV shows that the 
fan variables obtained with the least squares frictional loss factors gener- 
ally have smaller errors than those obtained with the point 13 loss factors. 

Table V lists the digital and hybrid. simulation steady-state data for point 
13 (Mach number - 0.8, 32 000 ft). For each simulation, data were taken with 
the frictional loss factors set to the design point values and with the fric- 
tional loss factors set to the values determined by the least squares fit. 

For the data presented, use of constant frictional loss factors (design 
point values or least squares values) give acceptable steady-state accuracy. 

If reasonable steady-state accuracy over the operating envelope Is desired, the 
least squares frictional loss factors can be used. However, If large pertur- 
bations from the operating conditions occur and/or good steady state accuracy 
Is required, then It would be necessary to Implement some functional form of 
the "exact" frictional loss factors. This would Increase the model complexity 
but should not Increase significantly the simulation computation time. 


Frequency Response Results 

For dynamic validation, a comparison with frequency response results from a 
distributed model of the AWT tunnel to the lumped parameter mathematical model 
was done. A distributed model of the AWT duct circuit was built using a method 
suggested In reference 2 In which subsonic duct sections are represented by 
linearized one-dlmenslonal wave equations. When applied to supersonic Inlets, 
this method was found to give good agreement with experimental data (ref. 3) 
and with a one-dlmenslonal method of characteristics approach (ref. 4). The 
equations relating conditions across the fan and the fan map were the same for 
both models. The distributed model was programmed on the same hybrid computer 
as the lumped parameter model. It was decided that the means of comparison 
would be a frequency response of test section Mach number to fan speed. 

Results were obtained at Mach 0.3 and 0.8 tunnel operating conditions (points 
3 and 13, table III) with only the pressure control loop closed. As mentioned 
earlier, the pressure control loop was closed to prevent the slow drift tend- 
ency of the hybrid simulation. Closing this control loop did not appear to 
affect the frequency response results. Also, the heat removed by the heat 
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fxchanger, Q^x. was Instantaneously set equal to the fan heat of compression, 
QpAN* Results for the Mach 0.8 condition are shown In figure 9 where normal- 
ized amplitude (ratio of change In Mach number to change In fan rpm, divided by 
the ratio at 0.05 Hz) and phase angle are plotted as a function of frequency In 
figures 9(a) and 9(b) respectively. The data were obtained by means of a sweep 
frequency technique. The disturbance variable, fan speed, was kept at a con- 
stant amplitude (10 rpm peak-to-peak) and the disturbance frequency was swept 
continuously from 0.05 to 10 Hz. Data beyond 4 Hz were discarded for reasons 
to be discussed In the next paragraph. The results show that the lumped param 
eter results have a greater amplitude for frequencies less than about 2.5 Hz 
and more attenuation for frequencies greater than 2.5 Hz. The phase data agree 
very well over the range of- 0.4 to 3 Hz. The 180° phase shift point, which Is 
Important from a controls pol nt-of-vlew. Is predicted to occur at slightly less 
than 2 Hz by both models. Also, the predicted amplitude ratios compare favor- 
ably In the area of 2 Hz (180* phase shift point). Similar agreement was 
obtained for the Mach 0.3 operating condition. 

In reference 5, an approximate relationship Is given for determining the 
number of model lumps required to satisfy a specified accuracy requirement. 
Rewriting this relationship, a maximum error can be determined for a model with 
a specific number of lumps at a frequency of Interest. Given the relationship 



where 

N represents the number of model lumps, 
f the frequency of Interest, 

4 the maximum error . 

Substituting for f 0 , the expression a/1, where a Is the speed of sound 
In air and 1 Is the length of the tunnel In feet, and further substituting 
for a, the expression. 
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where T Is taken to be the static temperature of air In the tunnel In °R, the 
relationship can be rewritten as 
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Using the frequencies mentioned above ( 1 . e . , frequencies of 4 Hz and 2 Hz), 
the error for the 4 Hz frequency Is 0.914 and the error for the 2 Hz frequency 
Is 0.228. This seems consistent with the results observed In figure 9. It was 
felt that the five lump model would be adequate for purposes of Investigating 
wind tunnel controls, and that Its use would likely result In a somewhat con- 
servative control design. 



CONTROLS APPLICATION RESULT* - 

Using the lumped parameter model of the circuit and the models he tun- 
nel controls and subsystems that were described earlier, transients ..re run to 
demonstrate that tunnel conditions could be changed In a stable manner with a'l 
control loops closed. Th ® tunnel design condition of Mach 0.8, 32 000 ft alti- 
tude was chosen as the operating point. The frictional loss factors were set 
to the design point values. Step changes of 10 percent In the total tempera- 
ture, total pressure, and Mach number setpoints were then made Individually to 
Investigate the controlled tunnel responses. The results obtained from the 
hybrid simulation, are shown In figures 10(a) to (c) for the total temperature, 
total pressure, and Hach number changes , respectively. Corresponding results 
from the digital simulation are shown In figure 11. 

Ihe setpoint change In test section total temperature. T , (f'gs. 10(a) 

TS 

and 11(a)), calls for an Increase from 456 to 501 C R. The following general 
observations are^made from the hybrid simulation traces. As expected, the 
amount of heat, Q^x. removed by the refrigeration system Initially goes down 
to allow temperature to Increase. The Increase ‘ ' f '-'oerature would also be 
expected to cause a decrease In Mach number a - ease In pressure, with 

out action by the Mach number and pressure co ,ps. As can be seen from 

the traces, fan speed Nf^ N did Increase to ; n Mach number and air was 

exhausted from the tunnel to maintain pressure, fomentary deviations can be 
seen for Mach number from 0.76 to 0.735 and for pressure (less than a 40 psf 
Increase). The tunnel Is seen to arrive at a new steady- itate condition within 
3 min with the fan operating at a slightly higher speed and the refrigeration 
system removing slightly less heat. A comparison of the heat removed transient, 
0 HX , shows a discrepancy between the hybrid (fig. 10(a)) and the digital 
(fig. 11(a)) results. The digital simulation result Is seen to settle out at 
a value somewhat higher than the predisturbance v.,1ue. This discrepancy 
between the hybrid and the digital results Is unexplained. However, as was 
pointed out earlier, a very small temperature difference in the models can 
make a significant difference In the fan heat of addition and hence, heat 
removed by the refrigeration system. It Is believed that the approximations 
In the hybrid model and/or analog component Inaccuracies could be the cause. 
Truncation and other errors associated with the digital simulation Integration 
algorithm could also be a contributing factor. Despite the discrepancy, no 
adverse control actions or setpoint Interactions are apparent during the 
transient. 

Ihe setpoint change In the test section total pressure (figs. 10(b) and 
11(b)), P-. , calls for an Increase from 870 to 960 psf. In this case supply 
TS 

air, m$A, Is Injected to raise the pressure. The traces show the alternate 
supply and exhaust air actions that are required to achieve the desired change 
In pressure. Control action by the temperature control loop Is also required. 
The amount of heat removed by the refrigeration system tan be seen to Increase 
to compensate for the additional air that Is being circulated In the tunnel. 

As expected, there Is very little Interaction with the Mach number. loop. This 
Is apparent from the sl’ght. momentary deviation In fan speed. The test sec- 
tion variables come to the new steady-state condition In approximately 2 min. 

The control action Is well behaved. 



The setpoint charge In test section Mach number (figs. 10(c) and 11(c)), 
*TS- calls for a decrease from 0.76 to 0.68. The main control variable for 
this transient Is fan speed n^n which Is seen to decrease from about 370 to 
355 rpm. There Is a corresponding decrease In heat removal Q^x by the 
refrigeration system and a slight addition of supply air m$A- The tunnel 
comes to equilibrium In about 1.5 min. 

Tne results of figures 10 and 11 Indicate that the tunnel controls are well 
behaved and that there are no control loop Interaction problems. Limited 
results at other tunnel conditions Indicated that variable controller gains 
and/or setpoint scheduling would be required to provide smooth transient oper- 
ation over the full operating range of the tunnel. The simulation, particu- 
larly on the hybrid computer, provides a useful tool for determining such 
schedules and for Investigating other tunnel control problems/sltuaticns . 


CONCLUDING REMARKS 

Computer simulations provide an Inexpensive and safe means for designing 
and predicting the performance of large test facilities. Simulations can be 
used to Investigate/assess alternative control schemes and operating scenarios 
to assure that the facility will operate safely, stably, and as Intended. 
Simulations also can provide a training ground for the facility operators. 

The lumped parameter wind tunnel model that was presented In this report. 

Is structured to allow the Incorporation of subsystems that provide functions 
ruch as refrigeration and supply/exhaust air. The model can, with some modi- 
fication, be used to represent a variety of closed-loop subsonic wind tunnel 
configurations. The model provides a useful analytical tool for predicting 
tunnel dynamic behavior and controls Investigation. It can also be made to 
accurately represent steady-state tunnel operation by means of friction loss 
factor terms Included In the momentum equations. The model can be simulated 
on either a digital or an analog computer. 

The model was used to simulate the proposed Lewis AWT. A comparison of 
steady-state results from a contracts supplied steady-state performance deck 
showed good agreement. A set of constant friction loss factors, determined by 
a least squares method, gave maximum errors of about 4 percent In total temper- 
ature. total pressure., and mass flow rate throughout the tunnel and over the 
tunnel's operating envelope. Generally the errors were about 1 percent. Large 
errors (on the order of 20 percent) In fan heat of compression and torque were 
observed. This Is due to extreme sensitivity to the temperature rise across 
the fan. A comparison of the tunnel's dynamic response predicted by the lumped 
parameter model to that predicted by a distributed (wave equation) model showed 
reasonable agreement. The conclusion was that the use of the lumped parameter 
model would result In conservatively designed controls. Transient results for 
the AWT model. Including transfer function models of the tunnel's refrigeration 
and supply/exhaust air subsystems and controls were presented. The results 
showed that, for step changes In test-section Mach number, total temperature, 
and total pressure, the tunnel response was stable and well behaved for the 
chosen control architecture. 



APPENDIX A 


SYMBOL TABLE 

A_ area, ft 2 

A_-_ average area, ft 2 

c v specific heat at constant volume, 0.17134 Btu (Ibm ®R) 

c p specific heat at constant pressure, 0.23988 Btu (lbm ®R) 

D_-_ average diameter, ft 

7_-_ frictional loss factor 

g c gravitational constant, 32.174 (lbm ft) (lbf sec 2 ) 

J mechanical equivalent of heat, 778.26 (lbf ft) Btu 

L_-_ length (AWT circuit lump), ft 

L_ length (specific section In contractor's steady-state performance deck). 

In. 

M_ Mach number 

m_ airflow, Ibm/sec 

Npan rotational speed of fan, rpm 

P_ pressure, lbf/ft 2 i 

Pj_ total pressure, lbf/ft 2 

Opan heat addition by fan, Btu/sec 

heat removal by refrigerator, Btu/sec 

1 

R gas constant of air, 53.34 (lbf ft) (lbm # R) j 

I 

T_ temperature, ®R j 

Tj_ total temperature, ®R 

i 

V_-_ volume (AWT circuit), ft 3 

V volume (specific section In contractor's steady-state performance deck), 

ft 3 

v_ velocity, ft/sec 

t Internal energy, Btu/ft 3 

y specific heat ratio, 1 .4 

14 

I 

l 


A 



n efficiency 

p density, lbm/ft 3 

t torque, 1 b f ft 

Superscripts : 

• time derivative 
average 
Subscripts : 

volume (station/section designation) 
lump 

a air 

c corrected 


EX 

exhaust air 

FAN 

fan 

HX 

heat exchanger 

PESA 

plenum air out 

PES8 

plenum air In 

SA 

supply air 
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APPENDIX 8 


EQUATIONS FOR CALCULATING TUNNEL VOLUMES, AVERAGE AREAS, 

LENGTHS ANO DIAMETERS 

To account for changes In tunnel diameter In a given lump (refer to 
fig. 2), average areas, average diameters, lengths, and volumes were precal- 
culated. Values for Inlet area, exit area, and length of each specific AWT 
section are shown In figure 81. Table BI contains the first part of the output 
from a database generation program. The following subsections give the 
equatlons/calculatlons for the areas, lengths, diameters, and volumes. 


VOLUME 


The sections of the AWT tunnel are considered to be frustoms of cones 
since the diameter of the tunnel varies throughout the section. Therefore, the 
volume for each section In figure Bl Is calculated from 


L_ 

V_ = — (A 1 ♦ A 2 ) 

(which Is equivalent to the equation for the frustum of a cone) 
where 


(B-l) 


V_ Is the volume calculated (see table Bl) 

L_ Is the length, L, given for each section In figure 81 

A^.Aj are the Inlet and exit areas respectively for each section In figure Bl 

The following equations show the calculations to obtain the six flow 
volume values for the AWT application. 


V, „ V, » V„ 

— - = 13 948 ft 2 


V v ♦ V ♦ V 

^ , - 3 2 4 5 - - 69 761 ft 2 


V 6-TS V 6 * V 7 * V 8 f V 9 * V 10 * V 11 * V 12 * 2 * V 13 


228 626 ft' 


V ip . 2 .11 = 25 4fc8 ft 2 


V v ♦ V ♦ V 

10 . 3 0 525 ft 2 


( B - 2 ) 
( B-3) 

( B- 4 ) 

(B-5) 

( 8 - 6 ) 
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Then 


1-3 2 

» 13 9<8 ft' 


- 1 -' 3 . * 83 309 ft 2 


8 2 


V , '-m . 398 387 ft 2 


3 2 


^5.^1 , 254 094 f. 2 


^06 . ^ ■ 55 ,93 ft 2 


■ —I -- = 30 525 ft 2 


LENGTH 

To calculate the length per lump (refer to fig. 2), the following equa- 
tions were used: 


L l-3 = h * L 2 * * + 3 - 15 = 35 • 5 ft 


L 3 _ 6 = L 3 ♦ l 4 ♦ L 5 = 24.5 ♦ 42.7 ♦ 0 » 67.2 ft 


(8-7) 

( 8 - 8 ) 


L 6-TS = L 6 + L 7 * L 8 * L 9 f L 10 * L 11 * L 12 * 2 L 13 


. 0 * 36.5 ♦ 70 ♦ 56.7 *• 0 ♦ 13.3 ♦ 47.7 ♦ £ (55) * 251.7 ft (8-9) 


1 


1 


■TS-16 = f L 13 * l 14 f L 15 = 2 (55) * 61 - 3 * 31 - 9 = 1?0 - 7 (B-10) 

L 16-19 = L 16 * L i 7 * L 10 = 36 » 38 > 15.7 = 87,7 ft (8-11) 


where l_ Is the length, L. of each section as shown In figure Bl. 


AVERAGE AREA 


Since the tunnel diameter changes through a lump (refer to figure Bl ) , an 
average area Is calculated for use In the momentum equation. 
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(B-12) 


•1-3 


1-3 * .2 ,2 

n 

V 1 V 2 


778.88 ft' 


■3-6 


3 ' 6 * l 2 l 2 l 2 

L 3 U 4 l 5 

V + V + V 
v 3 v 4 v 5 


1830.46 ft 


-6-TS 


6-TS 


L 2 L 2 L 2 L 2 L 2 L 2 L 2 L 2 ! 

L 6 L 7 L 8 U 9 L 10 L 11 L 1 2 \2 U 13/ 

— - * ♦ - — ♦ * -* — 


= 1186.37 ft' 


V 6 V 7 V 8 V 9 V 10 V 11 V 12 


- y 
2 V 1 3 


TS-1 6 


TS-16 n 2 ) 2 2 

\2 L 13/ i 4 L 15 

1 * y * V 

i V,, 'l 4 V 15 


= 411.08 ft' 


2 13 


'16-19 


16 * 19 * L ?6 . L ?7 . L ?8 
V 


693.05 ft' 


* y " f V 
16 17 v 18 


AVERAGE 0IAMETER 

Using the average area and the equation 

*• “*2 
4D‘ 

the average diameter can be calculated also for use In the momentum 
So rearranging 

7 * 

▼ 


0^3 = 31.491 ft 
0. , = 48.276 ft 

J-0 


(B-13) 

(B-14) 

(8-15) 

(8-16) 

( 8 - 1 7 ) 
equation. 

(B-18) 
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0 £ TC . 38.866 ft 

U“ \ O 

°TS-16 = ^ISJl 

0 16 . 19 - njoLU 
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TABLE BI. - TUNNEL GEOMETRY - DATABASE OUTPUT 


Section 

Inlet 

area 

Exit 

area 

Length 

Volume 


Fan exit straight 

600.84 

602.19 

4 

2406.1 

v ! 

Fan tal leone 

602.19 

1035.63 

31.5 

25 489.2 

V 2 

Tailcone strut 

1035.63 

1414.31 

24.5 

29 891.6 

V 3 

Cooler diffuser 

1414.31 

3930.37 

42.7 

109 630.7 

v 4 

Cooler 

3930.37 

3930.37 

0 

0.0 

v 5 

SI lencer 

3930.37 

4242.91 

0 

0.0 

v 6 

Cooler contraction 

4242.91 

2042.82 

36.5 

112 295.8 

v 7 

Corner 3 

2042.82 

2042.82 

56.7 

115 827.8 


v„ 

Crossleg 

2042.82 

2042.82 

13.3 

27 169.5 


8 

Corner 4 

2042.82 

2042.82 

56.7 

115 827.8 

V 9 

Honeycomb 

2042.82 

2042.82 

0 

0.0 

v io 

Stilling chamber 

2042.82 

2042.82 

13.3 

27 169.5 

V 11 

Test section contr. 

2042.82 

314.15 

47.7 

50 213.2 

V 12 

Test section front 

314.15 

322.07 

27.5 

8 747.7/ 


Test section back 

322.07 

329.99 

27.5 

8 965.51 

V 13 

High speed diffuser 

329.99 

506.23 

61 .3 

25 433.1 

V 14 

Corner 1 

506.23 

530.45 

31 .9 

16 533.5 

V 15 

Crossleg 

530.45 

759.10 

36 

23 089.2 

V 16 

Corner 2 

759.10 

710.99 

36 

26 456.9 

V 17 

Fan entrance duct 

710.99 

754.64 

15.7 

11 503.5 

V . Q 
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Total tunnel 



562.8 

736 656 

; 



Volume, 

vi-ni /2 

Length, 

L1-U1 

Average 
area , 

A1-W1 

Model volume 

13 948 

35.5 

778.88 

Model volume 

83 709 

67.2 

1830.46 

Model volume 

298 387 

251.7 

1186.3/ 

Model volume 

254 094 

120.7 

411.08 

Model volume 

55 993 

87.7 

693.05 

Model volume 

30 525 


_ . _ _ 

Total tunnel 

736 656 

562.8 




i 
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APPENOIX C 

CIRCUIT EQUATIONS AND CONSTANTS 


"'1-3 

*' " Vi 


(C-l) 


£ 1 V 1 


p, 2Jg c 

* - *’ 

’ ‘V^vT 

t t, ■ ! 1 (' * °- J "?) 


p r 'iVl 


p . p, (1 . 0.2 «;) 


Y/(Y-1 ) 





( C-2) 

( C-3) 

(C-4) 

(C-5) 

EC-6) 
( C-7) 

(C-8) 


A 1-3 9 c 


r3 


-1-3 


m l -3 

< P 1 * P 3> * x < V 1 - V 
9 c A l-3 


1-3 i -3 nil 
»e S l-3 ! 


(C-9) 


"'3-6 
3 “ '3*3 


C-l 0) 


!l V J- 
p 3 ■ 2Jg c 


(C-l 1 ) 
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-1 




T . T, (l . 0.2 H|) 


(C-12) 


( C-l 3) 


P 3 3 p 3 R a T 3 


■l. • "3 ( 


1 ♦ 0.2 HI 


V Y/(Y-1) 


p 3 ' jV } _ 2 V 3 _ 6 \ (m l-3 * m SA * m PES8 ' m 3-6 ) 


c 3 3 /V 13 P V 3 _A (m l-3 * m SA * m PESB )T 7 1 ' m 3-6 T T 3 

— - ♦ - — ) L 


(P 3- ‘v f rf” (v 3-v 2 

9 c A 3-6 g c°3-6 


6 ' p 6 A 6 


c Jb _ V 6_ 
p 6 ?Jg c 


K R a T 6 


\ . T^I * 0.2 H 6 
o 

V VA 


\ ■ p 6 0 • °- 2 "D 


2xy /( t-1) 
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1 


(C- 25 ) 


m 6 -TS = 


" ' V 3-6 VTs\ ( ™ 3 - 6 VTS) 


£fe | / V 3-6 VlS l 

V 2 '2 


m* / Ty m# Ty — 

3-6 T, 6 -TS T, c 
j bp 


HX' 


(C- 26 ) 


A 6 -TS q c 
L 6 -TS 


( p 6 - p ts ) . ^IS- ( , t - ,„) - 'to 


g c A 6 -TS 


6 TS' 


(C- 27 ) 


g c° 6 -TS 


m 


. TS— 1 6 
TS p TS A TS 


(C- 28 ) 


Its Its. 
p ts ' 2 Jg c 


'TS 


(C- 29 ) 


'TS 


TS 


V^cVtS 


T I ts • t ts (’ • 02 "?s) 


P TS = p TS R d T TS 


( C - 30 ) 

(C- 31 ) 
(C- 32 ) 


P I, S ■ P Ts(' * 0 2 "n) 


T/(T- 1 ) 


1 


PjS 1 V 6 -TS V TS -1 6 
2 * 2 


(m 6 -TS " m TS -16 ‘ m PESA^ 


( C- 33 ) 
(C- 34 ) 


TS I 6 -TS f V TS -16 


m 6 -TS T T 6 ‘ m PESA T T TS ‘ m TS- 16 T T 


J 


(C- 35 ) 
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TS-16 „ 
"TS-16 ’ L ts .„ 3 c 


>p .pi. „ , , . f TS-16 L TS-16 ^TS 

( TS 16 „ I TS V 1 6 r P TS 2 

9 c A TS- 16 g c°TS-16 

(C-36) 


"*16-19 


Vl6 ' P 16 A 16 


(C-37) 


' c 16 V 1 6 


p 16 2Jg c 


'16 


(C-38) 


'16 


1 6 




Vl6 


T l 16 = T 1 6 l 1 * °' 2 "? 6 ) 


P 16 1 p 16 R a T 16 


P I 16 = P 16 0 * °- 2 M h) 


t/(y-1) 


16 = /V Tr V, t ln \ {m TS-16 • m l 6-19 - "W 


'TS-16 16-19 

2 * 2 


( C - 3 9 ) 

( C - 40 ) 
(C-41) 

( C-42) 

( C-43) 


C 1 6 = /V 


»TS-16 ° V lb- I9 \ (*1S-16 T Tt s ' "EX'I, 6 ' *1 1 9 ' > , 6 ) 


.. A 1 6-1 9 

mi6 - 19 ’ L 16-19 


( p p ) , -16-19 , y i f 16 = .l_ 9 L 16-19 P 1 6 V 1 6 

' K 1 6 19 ; „ t * V 1 6 V 19 ; - 

9 C A 16-19 9 C°16-19 


2 

( C - 4 5 ) 


FAN 

3 19 A 19 


(C-46) 
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In hi. 

*19 " 2Jg c 


19 


(C-47) 


19 


19 


T T 1 T 19 
19 


V Y9 c R a T 19 

(’ * 02 "? 9 ) 


P 19 a p 19 R a T 19 


/ ? >y/(y-1) 

P T„ • P 19 (' * °- P "?9> 


( C - 48 ) 

( C-4S) 
( C-50) 

(C-51) 


*19 ’ /V 


1 


16-19 


(m 


16-19 " m FAN ) 


(C-52) 


19 /V 




(C-53) 


19 


'1 g * 2116.2 


(C-54) 


19 


T 19 518.69 


(C-55) 


FAN T 


19 


m FAN = f 


FAN 


FAN 


•v^ 


(C-56) 


C C - 5 7 ) 


"Van * "Van 


19 




(C-58) 


25 





T T * T 
TAN 1 9 



1 


1 .1 

/p\(lf-l)/Y 
P ' 1 -° 


n F AN 

v TAN 






Van * "Van ( t t fan - t t 19 ) c p 
x J 

t FAN = y FAN T 2^\~ 
y FAN 60 J 


T- 50) 
( C— 60 ) 
( C— 61 ) 
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APPENDIX C 


TABLE I. - AWT FAN MAP DATA 


0000630 


COnnON/TBLFl/NIK 3) ,N12C3) . 

N1313) ,N19(3),NF(3), ASTRQAt 52 ) , XM( 52 ) 

0000700 

c<*»» 

TABLE OF 

FAN MAP 

VALUES 





0 0 0 0800 

c * 

Ml = ARRAY OF MAP 

INFO FOR 

; MAP ROUTINES (I. 

E. MAP2 , 

MAPL) 

0 O 0 0 9 0 0 


DATA Ml/1 

,1,1. 7, 9/ 






0 0 1 0 0 0 

c» 

Y = FAN 

SPEED VALUES (9 CURVES) 




0 0 0 1 1 0 0 


DATA Y/100. ,150. ,225. ,250. ■ 

300. ,350 

.,900. ,997.67,500. 

/ 

000 1200 

c* 

X = FAN 

PRESSURE 

RATIO VALUES (7 

POINTS) 



0011320 


DATA X/ 

- 






0001920 


1 1 . 0005 . 

1.0030, 

1.0098, 

1.0066, 

1.0073, 

1.0096 , 

1.0108,- 

0 0 r, 1 9 9 0 


2 1.0017, 

1.0070. 

1.0130. 

1.0165, 

1.0201. 

1.0228, 

1.0296,- 

0-01960 


3 1 .0173, 

1.0332, 

1.0385, 

1.0951, 

1.0991, 

1 .0529, 

1.0596,- 

0001980 


9 1.0135, 

1 .0333, 

1 . 0989, 

1 . 0603, 

1.0696, 

1.0673, 

1 . 0686.- 

G-c : Sot 


5 1 .0519, 

1 . 0622, 

1 .0716, 

1.0858, 

1.0916, 

1.0963, 

1.1029.- 

0031510 


C 1.0790, 

1 . 0896 , 

1.0983, 

1.1175, 

1 . 1259, 

1.1319, 

1 .1903.- 

0 0 <i 1 59 0 


7 1.1003, 

1.1319. 

1.1597, 

1.1623, 

1.1721, 

1.1895, 

1.1890,- 

0001*60 


8 1 . 1299 , 

1 .1692 , 

1.1918, 

1.2186, 

1.2335, 

1.2393, 

1.2936,- 

•* 0 0 i 5 3 0 


9 1.2137, 

1.2525 . 

1.2793. 

1.2897, 

1.2990, 

1.3093, 

1.3127/ 

0 002 30 0 

' t 

Z = FAN 

CORRECTED 

FLOU VALUES CORRESPONDING 

TO FLOW 

VALUES 

tf'* 32 u 


DATA Z/ 

- 






1002090 


1 9986.6, 

9731.2, 

9975. 9, 

9258 ,9. 

3965.9, 

3980.5, 

2995.7,- 

1002360 


2 7511.7, 

7051.8. 

6502.7, 

6099.1, 

5999.2, 

5097.5, 

9986.1,- 

j.102380 


310519.5, 

9997.1, 

9113.9, 

8988.1, 

8028.9, 

7505.0. 

7083.9.- 

) Ou 2 CO 


912012.1, 

11019.8, 

9992.5, 

8983.5, 

8536.6, 

8000.6, 

7732.6,- 

n 0 C ? '• 2 0 


512989.2, 

12516 . 0 . 

12017.9, 

10995.3, 

10509.9, 

10011.9, 

9016.3,- 

7 0 o 2 9 9 0 


615052.9, 

19564 1, 

19029.0, 

13031.8, 

12533.9, 

12035.1, 

11026.5,- 

0002960 


717013.1, 

16027 .1, 

15059.8, 

19695.7, 

13981.2, 

13035.9, 

12550.6,- 

Of U?9S0 


819025. 1 , 

18037.9, 

17179.8, 

16066.8, 

15070.1, 

19585.7. 

19035.8,- 

s o o 


920059.3, 

19059.5, 

18086.8, 

17639.1, 

17089.5, 

16731.7, 

16080.2/ 


where 



FAN- 


FAN 


N 


FAN 


corrected fan airflow 

pressure ratio of fan 
corrected fan speed 


corresponds to 
z » f(x, y) 
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TABLE I. - FRICTIONAL LOSS FACTORS 


Design point 
values 

Least squares fit 
values 

^1-3 

0.02789 

f 1-3 

0.03017 

^3-6 

. 83294 

f 3-6 

.8402 

T 6-TS 

22.15813 

f 6-TS 

18.635 

f TS16 

.01574 

f TS-16 

.02076 

T 1 6-1 9 

.06969 

f 16-19 

.06858 
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TABLE II. - TUNNEL GEOMETRY CONSTANTS 


A. =. 600.34 ft 2 
V 

1-3 3 

-y 1 > 13948 ft J 

A - 778.80 ft 2 

*" 1 - 3 = 35.5 ft 

0. - = 31.491 ft 
i -3 

A 3 = 1035.63 ft 2 

, b37 09 ft 3 
A = 1830.46 ft 2 

L 3 _ 6 " 67 ' 2 ft 

0 6 * 48.276 ft 

A & = 3930.37 ft 2 

,398387 ft 3 

A 6 _ t = 1186.37 ft 2 
L 6-TS * 251 ' 7 ft 
0 6TS = 38.866 ft 
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TABLE II. - Continued 


A t$ * 322.07 ft 


a 254094 ft 3 


* TS -16 • 411.08 ft‘ 


L TS-16 * 120 - 7 ft 


°TS-16 - 22878 ft 


A 16 . 530.45 ft 


hsAk^h^l 3 55993 ft 3 


A lf> _ ig = 693.05 ft‘ 


*"16-19 = 67 ' 7 ft 


°1 6-19 - 29 ' 705 ft 


A ig = 754.64 ft‘ 


1 & : 1 - 9 - » 30525 ft 3 


31 



I 


TABLE III. - OPERATING POINT DEFINITION 
AT STATION 13 a 


Point 

Mach 

Altitude 

Total 

pressure 

Total 

temperature 

1 


762 

2116.798 

520.080 


m 

23 157 

875.587 

439.699 


.3000 


2116.799 

521 .798 




860.550 

439.670 


1 1 


2116.803 

524.184 


■ 

26 008 

840.147 

439.668 

7 


5 007 

2116.821 

527.905 

8 

-sMt '■ ' 9 

28 517 

809.517 

439.670 

mm 


14 875 

1531.682 

499.189 



30 478 

785.907 

439.668 

mm 


25 000 

1090.785 

471.724 

mm 


33 209 

753.555 

439.669 

Bn 


32 000 

875.939 

456.530 



55 000 

292.376 

439.886 

II 

.8367 

34 000 

828.063 

453.286 

mm 

.8367 

55 000 

303.381 

444.571 

19 

.9160 

54 999 

330.093 

455.411 


a Stat1on 13 Is entrance to test section. See 
figure B1 . 
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TABLE IV. - COMPARISON OF STEAOY-STATE ERRORS (PERCENTAGE ERRORS) 


FOR DESIGN-POINT AND LEAST-SQUARES FRICTION LOSS FACTORS 



Rotational 
speed of 
fan, 

n FAN* 

rpn 

Pressure 

ratio, 

( P/P) FAN 

Heat addition 
of fan, 

Qfan. 

8tu/sec 

Torque 
of fan, 

T FAN > 
lbf ft 

For all 
m. P T , Tj 


Point 3: Mach 

number, 0.3; altitude, 1710 ft 


Design 

Digital 

9.5 

1 .3 

48.2 

35.4 

1 .7 


Hybrid 

11.0 

1.0 

47.0 

27.0 

<2.0 

Least 

Digital 

4.4 

.7 

25.2 

20.2 

.9 

squares 

Hybrid 

6.4 

.2 

13.3 

5.3 

(a) 


Point 4: Mach 

number. 0.3 

altitude, 24 361 ft 


Design 

Digital 

8.4 

1.2 

41.3 

30.3 

1.7 


Hybrid 

11.0 

1.5 

42.0 

27.0 

<2.0 

Least 

Digital 

4.0 

.6 

20.9 

16.4 

.9 

squares 

Hyfcrld 

5.9 

.6 

17.6 

9.5 

(a) 


Point 8: Mach number, 0.52 

; altitude, 28 517 ft 


Design 

Digital 

4.5 

1.7 

24.5 

19.2 

2.6 


Hybrid 

6.0 

1.7 

24.0 

17.0 

<3.0 

Least 

Digital 

1 .2 

.5 

6.7 

6.0 

.8 

squares 

Hybrid 

1.4 

.3 

4.2 

2.6 

(a) 


Point 9: Mach 

number, 0.6 

altitude, 14 875 ft 


Design 

Digital 

1 .5 

1 .4 

14.0 

12.4 

2.4 


Hybrid 

4.4 

<2.0 

22.0 

17.0 

<3.0 

Least 

Digital 

,9 

.4 

5.0 

4.0 

(a) 

squares 

Hybrid 

1 .0 

.4 

5.0 

3.5 

(a) 


Point 13: Mach 

number, 0.8 

; altitude. 32 000 ft 


Design 

Digital 

(a) 

(a) 

(a) 

(a) 

(a) 


Hybrid 

(a) 

(a) 

(a) 

(a) 

(a) 

Least 

Digital 

-2.0 

-1.6 

-13.6 

-11.8 

-2.5 

squares 

Hybrid 

-2.2 

-1.7 

-14.4 

-12.7 

<3.0 


Point 16: Mach number, 0.8367; altitude, 55 

000 ft 


Design 

Olgltal 

-2.6 

-1 .6 

-11.8 

-9.5 

-1.1 


Hybrid 

-2.4 

-1.5 

-11.0 

-8.0 

(a) 

Least 

Digital 

-4.2 

-3.1 

-23.8 

-20.5 

-3.7 

squares 

Hybrid 

-4.8 

-3.3 

-23.1 

-19.5 

<-4.0 


a Less than I percent. 
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TABLE V. - MACH NUMBER, 0.8; ALTITUOE. 32 000 FT 



Contractor 
supplied 
steady-state 
perf ormarce 

Design point frictional 
loss factors 

Least squares frictional 
loss factors 


Digital 

Hybrid 

Digital 

Hybrid 

m 

6595.72 

6613.73 

6627.20 

6614.58 

6623.54 

”3 

6595.72 

6613.84 

6608.89 

6614.63 

6617.43 


6595.72 

6613.82 

6617.43 

6614.61 

6617.43 

^TS 

6595.72 

6613.73 

6619.87 

6614.59 

6623.54 

• 

m l6 

6595.72 

6613.73 

6617.43 

6614.58 

6617.43 

m FAN 

6595.72 

6613.73 

6614.99 

6614.57 

6621 .09 

\ 

p > 3 

916.38 

917.98 

917.36 

895.21 

894.62 

913.54 

915.13 

916.14 

892.14 

893.40 

p >. 

891.26 

892.80 

893.86 

869.11 

869.90 | 

0 

\ S 

P T 

873.01 

873.01 

872.19 

873.00 

1 

872.19 

834.62 

836.19 

834.66 

828.63 

827.64 

' 1 6 






P T 

813.68 

815.23 

816.34 

807.67 

808.87 

*19 






\ 

T T, 

474.80 

474.14 

475.52 

471 .68 

472.67 

474.80 

474.14 

475.49 

471.68 

472.52 

J 

\ 

456.53 

455.92 

457.36 

455.92 

456.45 

T t 

456.53 

455.91 

456.52 

455.92 

456.52 

t ts 






T t 

456.53 

455.91 

456.52 

455.92 

456.41 

T 1 6 
T t 

456.53 

455.91 

456.59 

455.92 

456.59 

T 1 9 






«1 

.2996 

.2998 

.3005 

.3074 

a .3082 

m 3 

.1681 

.1682 

a . 1 683 

.1722 

a . 1 721 

«6 

.0438 

.0438 

a . 0439 

.0450 

a . 045 1 

m ts 

.7578 

.7566 

.7578 

.7568 

. 7574 

m 16 

.3768 

. 3768 

a . 3772 

.3810 

a . 381 0 

Ml 9 

.2600 

.2601 

a . 2595 

.2628 

a . 262 1 

P 1 

861 .11 

862.51 

862.89 

838.42 

839.23 

P 3 

895.63 

897.24 

898.44 

873.86 

875.40 

P 6 

890.06 

891 .60 

892.49 

867.88 

869.14 

P TS 

596.60 

597.41 

597.99 

597.24 

597.69 


*Not exp 1 Iclty available, calculated from other values. 

34 









TABLE V. - Continued 


J 

Contractor 

supplied 

steady-state 

performance 

Design point frictional 
loss factors 

Least squares frictional 
loss factors 

Digital 

Hybrid 

Digital 

Hybrid 

P 16 

756.72 

758.13 

758.97 

749.65 

750.88 

P 19 

776.30 

777.78 

779.72 

769.82 

772.10 

T 1 

466.45 

465.77 

d 467 . 08 

462.93 

a 463.86 

T 3 

472.16 

471.47 

a 472.81 

468.90 

3 4 69 .74 

T 6 

456.36 

455.75 

a 457. 18 

455.74 

a 4 5 6 . 26 

t ts 

409.43 

409.08 

a 409 . 49 

409.06 

a 409 .53 

T 1 6 

443.93 

443.32 

443.89 

443.05 

a 443 . 53 

T 19 

450.44 

449.83 

a 450.52 

449.71 

a 4 50 . 40 

P 1 

.0346 

.03472 

.03469 

.03395 

.03395 

P 3 

.0356 

.03568 

.03564 

.03494 

.03495 

p 6 

.0366 

.03668 

.03663 

.03570 

.03572 

P TS 

.0273 

.02738 

.02739 

.02737 

.02739 

p 16 

.0320 

.03206 

.03206 

.03172 

.031 76 

P 1 9 

.0323 

.03242 

.03246 

.03209 

.03215 

V 1 

317.21 

317.06 

318.70 

324.23 

325.47 

v 3 

179.09 

179.00 

179.69 

182.81 

183.04 

v 6 

45.90 

45.88 

46.01 

47.14 

47.19 

V TS 

751 .96 

750.04 

751.28 

750.31 

751 .65 

V 16 

389.13 

388.90 

390.14 

393.10 

393.77 

V 19 

270.54 

270.36 

270.08 

273.72 

272.64 

£ 1 

a 2.8348 

2.8403 

2.8406 

2.7645 

2.7649 

c 3 

a 2 .9028 

2.9050 

2.9077 

2.8304 

2.8339 

c 6 

a 2.8619 

2.865o 

2.8686 

2.7894 

2.7930 

C TS 

a 2 . 2234 

2.2266 

2.2279 

2.2262 

2.2278 

c 1 6 

2.5308 

2.5321 

2.5354 

2.5059 

2.5098 

c 1 9 

a 2 . 540 1 

2.5457 

2.5507 

2.5206 

2.5275 

(P/P) FAN 

1.1263 

1.1260 

1.1257 

1 .1084 

7.1079 

n fan 

373.9 

373.40 

371.95 

366.57 

365.84 

n fan 

398.5 

398.27 

296.33 

390.93 

389.83 

7 tfan 

a 474.81 a 

474.14 

474.83 

471 .68 

472.30 

m FAN 

c 

16094 

16101 

16052 

16253 

16206 


a Not explicit* available, calculated from other values. 
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